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Abstract 

The paper deals with a spatially homogeneous and anisotropic universe filled with perfect fluid and dark 
energy components. The two sources are assumed to interact minimally together with a special law of 
variation for the average Hubble's parameter in order to solve the Einstein's field equations. The law yields 
two explicit forms of the scale factor governing the Bianchi-V space-time and constant values of deceleration 
parameter. The role of dark energy with variable equation of state parameter has been studied in detail 
in the evolution of Bianchi-V universe. It has been found that dark energy dominates the Universe at the 
present epoch, which is consistent with the observations. The Universe achieves flatness after the dominance 
of dark energy. The physical behavior of the Universe has been discussed in detail. 
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1 Introduction 

Recent observations like type la supernova (SN la) [T] — [5] , CMB anisotropy [5] [7] , and large scale structure [5] 
strongly indicate that our Universe is spatially flat and there exists an exotic cosmic fluid called dark energy 
(DE) with negative pressure, which constitutes about 70 percent of the total energy of Universe. It is an irony 
of the nature and is a puzzling phenomenon that most abundant form of matter-energy in the Universe is most 
mysterious. Many cosmologists believe that the simplest candidate for the DE is the cosmological constant (A) 
or vacuum energy since it fits the observational data well. During the cosmological evolution, the A-term has 
the constant energy density and pressure p( de ) = —p( de \ However, one has the reason to dislike the cosmological 
constant since it always suffers from the theoretical problems such as the "fine-tuning" and "cosmic coincidence" 
puzzles [5]. That is why, the different forms of dynamically changing DE with an effective equation of state 
(EoS), w( de ) = p( de ) 

J pl de ) < —1/3, have been proposed in the literature. Other possible forms of DE include 
quintessence (u/ de ) > — 1)[TO], phantom (u/ de ) < -I) [TT] etc. While the possibility u/ de ) << -I is ruled out 
by current cosmological data from SNe la (Supernovae Legacy Survey, Gold sample of Hubble Space Telescope) 
[SlIU, CMBR (WMAP, BOOMERANG) [HI [14] and large scale structure (Sloan Digital Sky Survey) [15] data, 
the dynamically evolving DE crossing the phantom divide line (PDL) (ui^ de ^ = —1) is mildly favored. Some 
other limits obtained from the observational results coming from SN la data (Knop et al [16) ) collaborated 
with CMBR anisotropy and galaxy clustering statistics (Tegmark et al [T7]) are -1.67 < u/ de ) < -0.62 and 
-1.33 < w (de) < -0.79 respectively. 

Following Berman [18 and Kumar and Singh recently Singh et al. [2U] proposed a special law of variation 
for the average Hubble's parameter in Bianchi-V space-time, which yields a constant value of deceleration 
parameter. Such a law of variation for Hubble's parameter is not inconsistent with the observations and is also 
approximately valid for slowly time- varying DP models. The law provides explicit forms of scale factors governing 
the Bianchi-V Universe and facilitates to describe accelerating as well as decelerating modes of evolution of the 
Universe. Models with constant DP have been extensively studied in the literature in different contexts (see, 
Kumar and Singh 19 and references therein). Most of the models with constant DP have been studied by 
considering perfect fluid or ordinary matter in the Universe. But the ordinary matter is not enough to describe 
the dynamics of an accelerating Universe. This motivates the researchers to consider the models of the Universe 
filled with some exotic matter such as the DE along with the usual perfect fluid. Akarsu and Kilinc [2Tj — [2"5] 
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have investigated Bianchi-I and Bianchi-III DE models with constant DP. Yadav and Yadav [21] have studied the 
role of DE with variable EoS in Bianchi type-Ill Universe evolving with constant DP. Kumar [25, 2B] has studied 
some isotropic and anisotropic models of accelerating Universe with DE and constant DP. Recently, Yadav et al 
[2"7] has presented LRS Bianchi-V Universe with DE characterized by variable EoS assuming constant DP. 

In this paper, we have considered minimally interacting perfect fluid and DE energy components with constant 
DP within the framework of a Bianchi-V space-time in general relativity. The paper is organized as follows. In 
Section 2, the models and field equations have been presented. The Section 3 deals with the exact solutions of 
the field equations and physical behavior of the models. Finally, the results are discussed in section 4. 



2 Model and field equations 

The spatially homogeneous and anisotropic Bianchi-V space-time is described by the line element 

ds 2 = -dt 2 + A 2 dx 2 + e 2ax {B 2 dy 2 + C 2 dz 2 ) 7 (1) 

where A , B and C are the metric functions of cosmic time t and a is a constant. 

We define a = (ABC)* as the average scale factor of the space-time (JT]) so that the average Hubble's 
parameter reads as 

H= a ~, (2) 
a 

where an over dot denotes derivative with respect to the cosmic time t. 

The Einstein's field equations in case of a mixture of perfect fluid and DE components, in the units 8ttG = 
c = 1, read as 

1 

Rij ~ ~~ — -*ii> (3) 

where Ty = + T- ( / e) is the overall energy momentum tensor with Ty" and T^ de ^ as the energy momentum 
tensors of ordinary matter and DE, respectively. These are given by 

T ( ; n) 1 = diag[~p {m \p {m) lP {m \p {m) ] = dza.g[-l,w (m) ,cj (m) ,w (m V m) (4) 

and 

T ( / e) 1 = diag[-p (de \p (de \p( de \ P ide *>] = diag[-l,cj^ e \^ de \^]p^ (5) 

where and p^ are, respectively the energy density and pressure of the perfect fluid component or ordinary 
baryonic matter while u/ m ) = p\ m ) / p\ m > is its EoS parameter. Similarly, p( de ) and p( de ) are, respectively the 
energy density and pressure of the DE component while I pi de ) is the corresponding EoS parameter. 

In a comoving coordinate system, the field equations , for the Bianchi-I space-time fl} , in case of Q and 
d5j) , read as 

| + § + i§-5 = -- im) P im) - - (d£ V de) , (6) 

g + ^ + §^-^ = - (m V m) - - (de V de) , (7) 

i + ^ + A^~A^ = - Jm)p(m) - L " ide)p{d ^ (8) 



The energy conservation equation y = yields 



p (m) + 3(1 + u^)p^H + p [de) + 3(1 + u^)p^H = 0. (11) 
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3 Solution of Field Equations 

Integrating (fTU|) and absorbing the constant of integration in B or C, without loss of generality, we obtain 

A 2 = BC. (12) 

Subtracting ([6]) from {7}, ([6]) from (0) from and taking second integral of each, we get the following 
three relations respectively: 

= d\ exp (xi J a~ 3 dt^j , (13) 
^ = d 2 exp y aT 3 dt^ , (14) 

^ = d 3 exp ^x 3 J a~ 3 d?J , (15) 

where d\ , x\ , d% , X2 , d 3 and X3 are constants of integration. 

From equations (fil3 |) - ([T5"j) and (fT2|) . the metric functions can be explicitly written as 

A(t) = a, (16) 



B(t) = ma exp \ l j a 3 dt^j 
C{t)=m~ 1 aexp f—l J a~ 3 dt 



(17) 
(18) 



where 

s/TT , (^2 + X 3 ) 
m=y , d 2 d 3 , 1 = 

with 

d 2 =d^ 1 , X2 = -x\. (19) 

In order to solve the field equations completely, first we assume that the perfect fluid and DE components 
interact minimally. Therefore, the energy momentum tensors of the two sources may be conserved separately. 
The energy conservation equation y = 0, of the perfect fluid leads to 

p w + 3(l+w w )p w fl = 0, (20) 

whereas the energy conservation equation T^ e ' lJ = 0, of the DE component yields 

^ e '+3(l + w (de) )p (<ie) ff = 0. (21) 

Following Akarsu and Kilinc [3T], we assume that the EoS parameter of the perfect fluid to be a constant, 
that is, 

a/-) = ^ = const., (22) 

while a/ de ) has been allowed to be a function of time since the current cosmological data from SNIa, CMB and 
large scale structures mildly favor dynamically evolving DE crossing the PDL as discussed in Section 1. 
Now integration of equation (l20l) leads to 

p(-) = Co a-^ 1+uCm) \ (23) 

where Co is a positive constant of integration. 

Finally, we constrain the system of equations with a law of variation for the average Hubble parameter in 
Bianchi- V space-time proposed by Singh et al. [20] , which yields a constant value of DP. The law reads as 

H = Da~ n , (24) 

where D > and n > are constants. In the following subsections, we discuss the DE cosmology for n ^ and 
n = by using the law (|24l) . 
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3.1 DE Cosmology for n^O 

In this case, integration of (f24f leads to 



a(t) = (nDt + ci)~ , 



where ci is a constant of integration. 
Inserting (J25J) into ((T6j) - ((TH1) , we get 



A(t) = (nDt + d)- 
l 



B(t) = m(nDt + c{) ™ exp 
C(t) = m~ 1 (nDt + ci)» exp 



D(n - 3) 

-I 



D(n - 3) 



(nDt + ci)~^ 
(nDt + ci) 1 ^ 



(25) 

(26) 
(27) 

(28) 



provided n ^ 3. 

Therefore, the model (TTJ) becomes 



ds z = -dr + T~ dx z + m z e 



2 2ax+kT~ 



dy z + m e 



-2„2ax-kT~ 



dz' 



(29) 



where T = nDt + C\ and k — D ^_^ ■ 

The average Hubble's parameter (if), energy density (p^) of perfect fluid, DE density (p^ 6 **) and EoS 
parameter (u/ de )) of DE, for the model (|2"9")l are found to be 



H = D(nDt + ci) -1 , 

3(l + ^( m >) 



3(1 + 1. 



p (m > = c (nDt + ci)- 

p (de) = 2,D 2 (nDt + Cl y 2 - l 2 (nDt + ci)~» - 3a 2 (nDt + ci)~» - c (n£>t + ci) - " 

Wel _ (2n - 3)D 2 (nDt + c x )- 2 - l 2 (nDt + ci)~» + a 2 (nDt + ci)~» - CoW (ro) (nDi + Ci)" 

3D 2 (nDt + ci)- 2 - Z 2 (n£>t + ci)~» - 3a 2 (nDt + ci)~s - c (n£>t + c x ) 

The above solutions satisfy the equation (|2"9"|) identically, as expected. 
The spatial volume (V) and expansion scalar (8) of the model read as 

V = a 3 = (nDt + ci)^, 

= 3H = 3D(nDi + ci) _1 . 
The anisotropy parameter (A) and shear scalar (a) of the model are given by 



(30) 

(31) 
(32) 

(33) 
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A B 

A~B 



B_C 
B~C 



C^A 
C~A 



2l 2 
3D 2 



(nDt + a)' 



a 2 = -AH 2 = l 2 (nDt + Cl y 



The value of DP (q) is found to be 



q 



aa 

-Ta = n-1, 



(34) 
(35) 

(36) 
(37) 
(38) 



which is a constant. The sign of q indicates whether the model inflates or not. A positive sign of q, i.e., 
n > 1 corresponds to the standard decelerating model whereas the negative sign of q, i.e., < n < 1 indicates 
inflation. The expansion of the Universe at a constant rate corresponds to n = 1, i.e., q — 0. Also, recent 
observations of SN la [1] — [5] reveal that the present Universe is accelerating and value of DP lies somewhere in 
the range — 1 < q < 0. It follows that in the derived model, one can choose the values of DP consistent with the 
observations. 

We observe that at t = —c\/nD, the spatial volume vanishes while all other parameters diverge. Therefore, 
the model has a big bang singularity at t = —ci/nD, which can be shifted to t = by choosing c\ = 0. The 
singularity is point type as the directional scale factors A(t), B(t) and C(t) vanish at the initial moment. 
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Figure 2: Matter energy density vs. time with c = 1.5, c\ = 1, D = 2, n = 0.5, a/ 1 ™) = 0. 
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Figure 4: EoS parameter of DE vs. time with c = 1, c x = 0.05, D = 2, n = 0.5, w( ,n > = 0, / = 1.2, a = 0.7. 



1.2 
1 

0.8 
0.6 
0.4 
0.2 

0, 



- fiC*0 

— n 



o 



0.2 



0.4 



0.6 



0.8 



Figure 5: Density parameters vs. time with D = 2.5, cq = 12, ci = 1, u/ m ) = 0, n = 0.05, / = 1.2, a = 0.7. 



Fig.l shows that the directional scale factors monotonically increase with time. Thus, expansion of the 
Universe takes place in all the three spatial directions. 

From Fig. 2 and Fig. 3, we observe that p^> as well as p( de ) remain positive during the cosmic evolution. 
Therefore, the weak energy condition (WEC) as well as null energy condition (NEC) are obeyed in the derived 
model. Further, p( m ) and decrease with time, and approach to small positive values at the present epoch. 

The parameters H , 9 and a 2 start off with extremely large values, and continue to decrease with the expansion 
of the Universe. The anisotropy parameter A also decreases with the cosmic evolution provided n < 3. This 
shows that anisotropy of the model goes off during the cosmic evolution. The spatial volume grows with the 
cosmic time. 

Fig. 4 depicts the variation of EoS parameter (u/ de )) versus cosmic time for accelerating phase of Universe 
(q = —0.5), as a representative case with appropriate choice of constants of integration and other physical 
parameters . The SN la data suggests that —1.67 < Lo^ de ^ < —0.62 [TB] while the limit imposed on u/ de ) by a 
combination of SN la data (with CMB anisotropy) and galaxy clustering statistics is —1.33 < u/ de ' < —0.79 
[17]. Fig. 4, clearly shows that u! t - de ' evolves within a range, which is in nice agreement with SN la and CMB 
observations. 

The perfect fluid density parameter (fl^) and DE density parameter (fj( de )) are given by 

rt m} = ■ffxinDt + cx) 2 "'' 1 ^""'' , (39) 
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n(de) = 1 " ^ {nDt + cl) ^ - ^ {nDt + cl) ^ - ^ {nDt + C1 ) 2 "' <lr(m)) - ( 4 °) 

The overall density parameter (Q) is given by 

/2 2 

Q = fi< m > + r» (de ' = 1 - ^(nDt + d) 2 ^ 2 - ^(n,Dt + ci) 2 ^ 2 . (41) 

From equation (|41[) . it is observed that for < n < 1 , the overall density parameter (f2) approaches to 1 at 
late times. Thus, the derived model predicts a flat Universe at the present epoch. Fig. 5, depicts the variation 
of density parameters versus cosmic time during the evolution of Universe. It is observed that initially Universe 
was matter dominated and later on DE dominates the evolution of Universe which is probably responsible for 
the accelerated expansion of present-day Universe. 



3.2 DE Cosmology for n = 

In this case, integration of (1241 yields 

a{t) = c 2 e Dt , (42) 

where c 2 is a positive constant of integration. 
The metric functions, therefore, read as 

A(t) = c 2 e Dt , (43) 
B(t) = mc 2 exp [ot - ^L_ e ~ ZDt ^j , (44) 

C{t) = m- 1 c 2 exp ( Dt + ^r^e- 3Dt ) . (45) 
V 3Dc 2 ) 

Therefore, the model (JXJ) becomes 

ds 2 = -dt 2 + c 2 2 e 2Dt (dx 2 + m 2 e 2ax+hle ~ 3Dt dy 2 + m^e 2 ^-^ ' dz 2 ) . (46) 

where h =3lk- 

The average Hubble's parameter, energy density of perfect fluid, DE density and EoS parameter of DE, for 
the model (|46l) are obtained as 

H = D, (47) 
pim) = CoC - 3 ( 1 +- (m) ) e -3^( 1 +- (m) ) t , (48) 

p (de) = 3D 2 _ f c -6 e -6Dt _ io 2 c -2 e -2Dt _ ^-H ) e -3D(lW"" )t _ (4g) 

_on2 _ ;2 -6 -6Dt , n 2 -2 -2Dt 

{de) = 3D l_c^e +a c 2 e 



3L> 2 - l 2 c^ 6 e- 6Dt - 3a 2 c 2 2 e- 2Dt - c c 2 3{1+uj(m)) ' e -i">(iW»J)t ' 

The above solutions satisfy the equation (|4l)|) identically, as expected. 
The spatial volume and expansion scalar of the model read as 

V = 4e 3D \ (51) 

6 = 3D. (52) 

The DP is given by 

q=-l. (53) 
The perfect fluid density parameter and DE density parameter are given by 

-3(lW m) ) 

n (m) = £0£2 e -3D(l+ W M) tj (M) 



3L> 2 

;2 -6 2 -2 -3(lW m) ) 
o( de ) — 1 _ 1_£2_ -6Dt _ a c 2 „-2Dt _ C C 2 -3D(l+w (wi) )t 

~ 3D* 6 D 2 3D 2 6 



(55) 
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The overall density parameter reads as 



72 -6 2 -2 

Q = 1 -^ e ~ 6 -^H (56) 



The anisotropic parameter and shear scalar of model (1461) are given by 

o/2 -6 

a 2 = l 2 c- 6 e' 6Dt . (58) 

Recent observations of SN la [T] — [5], suggest that the Universe is accelerating in its present state of evolution. 
It is believed that the way Universe is accelerating presently; it will expand at the fastest possible rate in future 
and forever. For n = 0, we get q = — 1 ; incidentally this value of DP leads to dH/dt = 0, which implies the 
greatest value of Hubble's parameter and the fastest rate of expansion of the Universe. Therefore, the derived 
model can be utilized to describe the dynamics of the late time evolution of the actual Universe. So, in what 
follows, we emphasize upon the late time behavior of the derived model. 

From Equation (IBTfl) . we observe that o/ de ) « —1 for sufficiently large time t. Therefore, the late time dynam- 
ics of EoS parameter ui^ de ^ represents the vacuum fluid dominated Universe, which is mathematically equivalent 
to cosmological constant. Further, at late times, we have 

p( m ) w 0, 

p( de ) « 3D 2 , 

Asfl. 



Thus, the ordinary matter density becomes negligible whereas the accelerated expansion of Universe continues 
with non-zero and constant DE density, which is consistent with recent observations. Also the late time dynamics 
of the derived model shows that the anisotropy of the Universe damps out and Universe becomes isotropic. From 
equation (|56l) . it is observed that for sufficiently large time, the overall density parameter approaches to 1, i.e., 
f2 ps 1. Thus, the model predicts a flat Universe. 



4 Concluding Remarks 

In this paper, we have investigated the role of DE with variable EoS parameter in the evolution of Universe 
within the framework of a spatially homogeneous and anisotropic Bianchi-V space-time by taking into account 
the special law of variation of average Hubble parameter, which yields a constant value of DP. The main features 
of the work are as follows: 



• The models are based on exact solutions of the field equations. 

• In the present models, the matter energy density and DE density remain positive. Therefore, the weak and 
null energy conditions are satisfied, which in turn imply that the derived models are physically realistic. 

• The singular model (n ^ 0) seems to describe the dynamics of Universe from big bang to the present epoch 
while the non-singular model (n = 0) seems reasonable to project dynamics of future Universe. 

• The age of Universe, in the singular model, is given by 




which differs from the present estimate, i.e., To = H 1 ss \AGyr. But if we take D = 1, the model (n ^ 0) 
is in good agreement with the present age of Universe. 

• The EoS parameter of DE evolves within the range predicted by the observations. 

• The DE dominates the Universe at the present epoch, which may be attributed to the current accelerated 
expansion of the Universe. 
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• The Universe acquires flatness with the dominance of DE (see, Fig. 5). Thus, flatness of the Universe 
seems a natural consequence of DE. 

Hypothetical DE is the most popular way of explaining why the Universe is expanding at an ever increasing 
rate. DE plays a massive part in shaping our reality however nobody seems certain of what the dang stuff 
actually is. Future space missions hope to solve this mystery and shake up our current understanding of the 
Universe. 
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